Quantitative phase reconstruction for orthogonal-scanning 
differential phase-contrast optical coherence tomography 
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We present differential phase-contrast optical coherence tomography (DPC-OCT) with 
two transversally separated probing beams to sense phase gradients in various directions 
by employing a rotatable Wollaston prism. In combination with a two-dimensional mathe- 
matical reconstruction algorithm based on a regularized shape from shading (SfS) method 
accurate quantitative phase maps can be determined from a set of two orthogonal en-face 
DPC-OCT images, as exemplified on various technical samples. Original paper: Optics 
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Optical Coherence Tomography (OCT), a depth- 
resolved imaging method for the characterisation of 
turbid media, has become a well-established tech- 
nique in biomedical diagnostics and research. In case 
of transparent and weakly scattering specimens - like 
single cell layers - it was shown that enhanced con- 
trast can be obtained by combining the OCT prin- 
ciple with that of differential interference contrast 
(DIC) microscopy, resulting in differential phase con- 
trast OCT (DPC-OCT) [1]. With DPC-OCT, sub- 
wavelength differences in optical path lengths (OPL) 
could be quantitatively measured [2, 3], even below 
extended scattering layers. However, the phase dif- 
ference is determined between two adjacent probing 
beams, leading to unidirectional image contrast along 
the direction of beam separation and leaving features 
with phase gradients perpendicular to the separa- 
tion direction undetected. In addition, the unidirec- 
tional contrast in the images unavoidably causes se- 
vere streaking artifacts when quantitative phase maps 



are obtained by numerical line integration, as pro- 
posed in [3], due to noise accumulation and undeter- 
mined integration constants. These artefacts can be 
avoided if a well defined fiat reference plane is avail- 
able close to the structures of interest: a collinear, 
phase-sensitive approach in the time-domain (TD) 
was used to determine the phase difference between 
such an additional reference plane and the sample it- 
self [3]. However, this technique requires that even 
for the evaluation of a single sample layer a whole 3D 
volume - including both, layer and reference plane - 
has to be acquired. 

For the recently developed phase-sensitive spectral- 
domain (SD) OCT techniques a flat reference plane 
is equally placed in the vicinity of the sample struc- 
tures to form a stable common path interferome- 
ter [4, 5]. Although these SD-methods outperform 
their TD-counterparts [1, 2, 3] in terms of phase 
stability, TD-OCT techniques are still favourable in 
certain cases, like transversal scanning and full-field 
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(FF) OCT [6, 7] in combination with dynamic fo- 
cusing leading to a constant high lateral resolution 
throughout the whole imaging depth. In addition, 
phase retrieval with SD-techniques requires the acqui- 
sition, calculation and evaluation of whole 3D-data 
sets: three orders of magnitude more data is needed 
for SD-OCT compared to TD-OCT in order to ob- 
tain an en-face scan (C-scan) of a single interface. 
Consequently, acquiring individual C-scans is still 5- 
10 times faster with TD-techniques [6] despite the 
high acquisition rates of current SD methods [8]. 

We present a novel method which combines or- 
thogonal scanning DPC-OCT in the TD with a 
two-dimensional phase reconstruction approach and 
which is fully compatible with high-speed transversal 
scanning OCT as well as with FF-OCT. Our method 
effectively eliminates streaking artefacts so that ac- 
curate quantitative phase maps can be obtained, as 
exemplified on diverse technical samples, underlin- 
ing also the need of enhanced contrast and quanti- 
tative information retrieval for alternative OCT ap- 
plications outside the biomedical field [9]. 

The DPC-OCT setup is schematically depicted in 
Figure 1 and has been realised by modifying an exist- 
ing polarisation sensitive (PS)-OCT system for strain 
imaging in polymers and composite materials (cen- 
ter wavelength A=1550 nm, 90 dB system sensitiv- 
ity, 19 fim axial resolution in air) [10]: a Wollas- 
ton prism (WP1) is inserted in the sample arm in 
the focal point of the achromatic lens L with focal 
length / = 4.5 mm. WP1 has a divergence angle 
a = 5', thus a lateral shear r = 6.5 /urn, determining 
the lateral resolution, is introduced between the two 
polarisation-encoded probing beams (Figure 1 (b)), 
as verified by calibration measurements on a sharp 
reflecting edge. In addition, WP1 can be freely ro- 
tated around the beam axis, so that the shear can 
be oriented both in x or y direction as t x = Ax 
and T y = Ay, respectively. Due to refractive in- 
dex changes or thickness variations an OPL differ- 
ence causes a phase shift Aipi (with i = x, y) be- 
tween the two sheared beams. In the backreflected 
beam this generally results in an elliptical polariza- 
tion state which is brought to interference with the 
linearly polarised light of the reference beam (ori- 
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Figure 1: (a) Scheme of the DPC-OCT setup with 
indicated polarization directions, (b) configuration of 
the probing head with rotatable Wollaston prism (/ 
focal length, a divergence angle, Ax probing beam 
shear) 

ented at 45°). The two depth (z) resolved interfer- 
ence signals S\(x, y, z) and S%(x, y, z) are detected at 
each lateral position (x, y) in the PS-detection unit 
(WP2 and two photodetectors Di and D 2 ). Scaling 
the directional phase shift Aipi with the beam shear 
Ti yields the phase gradient tpi, 

V i { X ,v) = ^= miX ' V) - ,PliiX ' V \ i = { X ,v). (1) 

n ^ 

The two phase functions tpu and tf 2i of the detector 
signals Si(x,y,z) and z) are determined by 

applying the Hilbert transform H [2] with 



^>i,2x{x, y) = arctan 



H(S h2 (x±T x /2,y,z)) 
Si, 2 {x ±T x /2,y, z) 



(2) 



for the beam shear in x-direction and for tpi t 2 y in- 
troducing the shear in y-direction by an analogous 
expression. 

For the estimation of the achievable accuracy of 
the setup, we determined the phase distribution by 
repeating 2500 A-scans on a simple glass slide and by 
evaluating the 40 dB reflection of the bottom surface. 
Despite the rather long interferometer arms (>0.5 m), 
the standard deviation of the phase fluctuations was 
determined to be only 0.0055 rad corresponding to 
0.7 nm OPL at a center wavelength of 1550 nm. 

In Figure 2 Aipi(x,y) maps are depicted for a mi- 
crolens array (Suss MicroOptics, Switzerland) with 
known dimensions in order to validate our method 
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Figure 2: (a) Sketch of a microlens array test struc- 
ture with 4 interfaces, (b) cross-sectional DPC-OCT 
image, (c, d) en-face DPC-OCT images in x- and y- 
direction at interface (2), (e, f) en- face DPC-OCT 
images at interface (4), (g) reconstructed cumulative 
phase map at interface (2), (h) cumulative phase map 
at interface (4). White scale bar: 100 ^m 

(Figure 2 (c)-(f)). The array was placed upside down 
on a microscope slide to simulate several buried layers 
(cf. cross-sectional sketch and corresponding DPC- 
OCT scan in Figure 2 (a), (b)). Using both orthogo- 
nal DPC images Aip x (x, y) and Aip y (x, y) of the indi- 
vidual layers, the cumulative phase map <p{x, y) can 
be mathematically reconstructed with a regularized 
shape from shading (SfS) approach [11], which we 
adapted for the first time for the phase reconstruc- 
tion from DPC-OCT images (Figure 2(g),(h)). In 
this approach the energy functional £ according to 
equation (3) has to be minimized with respect to the 
total phase <p including a penalty term with the reg- 
ularization parameter \x for higher order derivatives, 

£{<t>{x,y)) = I {{f(x,y) -TL(0 x ,ip y )f 
Jn 

+Kf )2 xx+2 ( Ply + f )2 yy), (3) 

with the measured data f = I x I and R = 

\ fv J 

(^ x | denoting the derivatives of the total phase 
0y J 

(p to be reconstructed. In case of the microlens ar- 
ray the mathematical reconstruction of the OPL map 
at interface (2) yields - for a sample refractive index 
n s =1.449 at A=1550 nm - a total geometrical height 
of the individual lenslets of 1.64 fim, closely matching 



the technically specified value of 1.63 /im. 

In Figure 3 the performance of the SfS algorithm 
is compared to other reconstruction methods, using 
Aip x (x,y) and Aip y {x,y) maps of the top reflection 
of a thin corrosion inhibiting oil layer with small em- 
bedded abrasion particles. At first, we tested straight 
forward numerical line integration on a single en- 
face DPC image (Figure 3(b)), which leads to the 
above described streaking artifacts (Figure 3 (c)), es- 
pecially severe for samples with low signal-to-noise ra- 
tio (SNR). Using both DPC scans (Figure 3 (a), (b)) 
we then reconstructed the OPL map with a global 
Fast Fourier Transform (FFT)-based approach by 
solving the Laplace problem in the Fourier domain 
[12]. This method performs better than line inte- 
gration, especially for samples with low SNR (Fig- 
ure 3(d)). Nevertheless, slight blurring and haze is 
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Figure 3: DPC-OCT image of an oil layer with em- 
bedded particles, (a, b) DPC images with shear in x- 
and ^/-direction, respectively, (c) OPL reconstruction 
by numerical line integration using the single scan 
displayed in (b), (d) FFT-based OPL reconstruction 
from scans (a) and (b), (e, f) OPL reconstruction with 
SfS method using scans (a) and (b). White scale bar: 
100 fim 

present, partly caused by the low pass characteristics 
of the inverse Laplace filtering. Boundary artefacts, 
inherent to FFT-based methods, can be reduced by 
image continuation and windowing, but will at the 
same time deteriorate the results for extended sample 
features. In comparison, reconstruction by SfS yields 
a significantly improved result as depicted in Fig- 



3 



ure 3(e),(f) with no streaking artefacts or blurred 
features and areas. 

As a last example we recall that the measured 
phase differences are mapped into the interval [— n, ir] 
and that phase unwrapping is required in case of steep 
phase gradients as depicted in Figure 4 (a) , (b) for a 
soft and uneven acrylic polymer film. We applied a 
quality guided unwrapping approach [13] to obtain 
the two orthogonal unwrapped DPC en-face scans 
(Figure 4 (c), (d)). The SfS reconstruction of the un- 
wrapped images leads to the cumulative phase distri- 
bution (Figure 4 (e)), demonstrating that this method 
- also in combination with phase unwrapping - acts as 
robust method coping with steep and extended fea- 
tures as well as handling low SNR images like those 
in Figure 3. 




Figure 4: (a, b) Wrapped DPC images of a transpar- 
ent polymer layer (shear in x- and y-direction, re- 
spectively), (c, d) unwrapped phase of (a) and (b), 
respectively, (e) total phase map reconstructed from 
(c) and (d). White scale bar: 40 /xm 

In conclusion, we have demonstrated TD-DPC- 
OCT imaging of en-face phase gradients in variable 
directions combined with a regularized SfS recon- 
struction approach to determine quantitative OPL 
maps: from two orthogonal shear directions a more 
robust solution is obtained compared to conventional 
reconstruction methods, like numerical integration 
of single DPC-OCT images. Whereas the current 
setup has been realised with a rotatable Wollas- 
ton prism, polarisation-dependent electro-optic beam 
steering could be used in the future. Furthermore, our 



method is fully compatible with high-speed transver- 
sal scanning as well as with FF-OCT and provides en- 
hanced contrast and quantitative phase information 
for applications ranging from biomedicine to indus- 
trial metrology. 
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